Electron spin resonance shift in spin ladder compounds 
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We analyze the efTects of diflferent coupling anisotropies in spin-1/2 ladder on the Electron Spin 
Resonance (ESR) shift. Combining a perturbative expression in the anisotropies with temperature 
dependent Density Matrix Renormalization Group (T-DMRG) computation of the short range cor- 
relations, we provide the full temperature and magnetic field evolution of the ESR paramagnetic 
shift. We show that for well chosen parameters the ESR shift can be in principle used to extract 
quantitatively the anisotropies and, as an example, discuss the material BPCB. 
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Quantum magnetism is a fascinating branch of solid 
state physics [l[, and more recently of quantum gases in 
optical lattices [2[. Due to the competition between the 
exchange interactions between spins localized on different 
lattice sites, many fascinating quantum phases can occur 
ranging from ferromagnetically or antiferromagnetically 
ordered states to spin liquids. Recently spin systems 
have also allowed one to address questions usually re- 
lated to itinerant quantum systems such as Bose-Einstein 
condensation (EEC) [3| and Tomonaga-Luttinger liquid 
physics y. 

For potential application it is particularly important 
to have efficient probes to characterize the properties of 
magnetic systems. Fortunately powerful probes, such as 
neutron scattering [5| or nuclear magnetic resonance [6| 
exist and give important physical information. However 
in systems which are particularly sensitive to the sym- 
metry of exchange interactions (such for example the 
ones where BEG can be reahzed), a direct probe of small 
asymmetries of the coupling constants is desirable. An 
ideal probe for this purpose is the Electronic Spin Res- 
onance (ESR) which allows one to investigate the fre- 
quency shift of the paramagnetic resonance induced by 
small anisotropies |7H9|. In order to interpret the exper- 
imental observations an efficient theoretical analysis of 
the corresponding spectra is mandatory. Although the- 
oretical interpretations have been developed for a long 
time BiE 
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their range of applicability is generally 
limited to several models in a restricted domain of pa- 
rameters: low temperature region of spin-1/2 antiferro- 



magnetic chains 1^, high temperature limit of spin-1/2 



antiferromagnetic chains |13j . low temperature region of 
spin-1 Haldane chain compound [IJ, ll5|, spin-1/2 chains 
with nearest and next nearest neighbor interactions [16| , 
strongly dimerized spin-1/2 ladders analysed in a single 
rung picture [l7|. 

In this letter we present a flexible theory of the ESR 
spectra. Gombining an analytical expression of the ESR 
shift with a numerical evaluation of short range correla- 



tions we give the full dependence of the ESR shift with 
the magnetic field and the temperature. We apply this 
approach to the case of spin-1/2 ladders using temper- 
ature dependent density matrix renormalization group 
(T-DMRG) [11, [ii]. We show in particular that the ESR 
technique can be used to distinguish the nature of differ- 
ent anisotropies quantitatively. We illustrate our theory 
by applying it to ESR measurements on the compound 
(G5Hi2N)2GuBr4 (BPGB) that have been performed [13 
recently. Our study extends the analysis of the ESR 
measurements from Ref. [l7| focussing on the temper- 
ature and magnetic field dependence of the shift of the 
paramagnetic peak. We discuss the nature of different 
anisotropies deduced for BPCB and make proposals for 
additional experiments. 

A spin-1/2 ladder in a magnetic field with interaction 
exchange anisotropies can be described by the Hamilto- 
nian "H = THP -\- H' which consists of two parts. The first 
part 
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J\\ ^S'j,r5'i+i,; + Ji^S'i,i-5j,2-5M-B^'5'^ (1) 



includes the isotropic exchange interaction and the Zee- 
man term. The spin operator is defined by S',,; = 
{Sfi,Sfi,S^i) ~ <T, where a is the Pauli matrix act- 
ing on the spin located on leg ^ = 1,2 and rung i, and 
S = J2i I '^i,i is the total spin. Here Ji\ , J± > are the 
antiferromagnetic couplings along the legs and the rungs 
of the ladder, g is the Landc factor, fig is the Bohr mag- 
neton and h is the strength of the static magnetic field. 
Note that by this giigh has the unit of energy. 
The anisotropic exchange interaction 



"^ ^ Z^ 'h^P^i-l^i+lA "^ Z^ J±,p^i,l^i,2 



(2) 
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contains different spatial anisotropies J',, and J^ in the 
exchange interactions along the legs and the rungs, re- 
spectively, p = a,b,c represents the principal axis coordi- 
nates aajba, Ca of more general interactions A"^ S'^i^'^k 



instead of J'^^SfiS^ ^. in ^ (for dipolar interaction typ- 
ically) with A'^ , a =11, _L, arc symmetrical tensors in 
every coordinate basis x,y,z (denoted q,q' = x,y,z) 
aligned with the magnetic field h = hz. In the aa,ba, c^ 
coordinates, one denotes z = {za^a, Za^b, Za,c)- 

The ESR frequency shift 6uj is defined as the difference 
between the ESR paramagnetic frequency resonance, and 
the Zeeman frequency: 



hduj = hiUr — H. 



(3) 



Wr is the resonance frequency and H = 

Assuming that the anisotropics are small, i.e. 

I ph I>^1 pi '^ J\\:J±^ we determine ^ based on a first 
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order perturbation theory in the anisotropy. At this or- 
der, the ESR shift is represented by a static correlation 
function 041], 



hSio = — 



m\s+],s-])o 

2(5-)o 



(4) 



with S^ = 5*^ ± iS^ . The average {■)o is taken with 
respect to the unperturbed Hamiltonian Ti'^ ([1]). The 
formula Q holds for arbitrary temperatures T and mag- 
netic fields. For the spin ladder, the shift Q can be 
decomposed in two components 



hSu;^ J2 Uz)Y^{T,H) 
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corresponding to the effect induced by the anisotropy in 
the parallel (a =||) and the perpendicular {a ~±) in- 
teraction exchange. The temperature and magnetic field 
dependence of the shift ([5]) are determined by the short 
range correlations Ya (T, H) : 

Yii(r,g)= ^^^^^^^';^y^^^^-^^° , (6) 
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Computing these correlations is in general involved [20| . 
Here we use the T-DMRG method which allows us to 
determine the ESR shift. Note that this approach is quite 
general and can be extended to other systems. 

For the considered spin ladder we fix the direction of 
the magnetic field z with respect to the system orien- 
tation, the factors fa{z) depend only on the anisotropic 
coupling constants J^ : 



Uz)^t, 
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with t|| = 1 and t±^ = 1/2. Assuming that all the isotropic 
coupling is included in H°, i.e. J^p^'a.p = Oj ^^ ob- 




Figure 1. (color online): Magnetic field dependence of Y^ 
(lines with circles) for several temperatures and coupling ra- 
tios j. In (a-c), Y± is represented by lines without circles. 
In (d), the strong coupling j — ^ oo approximation ((9]) of Vn 
is shown with lines without circles. For the BPCB coupling 
ratio j — 3.6 in (c), the plotted temperatures are T ~ 0.72 K, 
T f» 1.44 K, T f» 3.6 K, T ^ 10.8 K in experimental unit. 
The two vertical dotted lines locate the critical fields Hd and 
Hc2 respectively. 



anisotropy A"^ along the magnetic field orientation z 
and has extrema when z \\ aa,ba or Cq. 

In the following we first discuss the properties of 
Ya{T,H) for different coupling ratios j = J\\/J±. We 
then focus on the compound BPCB and analyze its ESR 
measurements. 

At low magnetic field H ^ Hd both Y|| and Y± grow 
linearly with H (Fig. [J), i.e. Y\i,Yj_ (x H {Hd < Hc2 
are the two critical fields between which the system is 
gapless [19[). In a strong magnetic field H :$> Hc2 the 
system becomes totally polarized and y||,KL ~ 1/2. Be- 
tween these two limits the behavior of Yji and Y± is very 
subtle and clearly depends on the coupling ratio j. In 
the zero interladder coupling limit j = 0, the two ladders 
decouple. Y± reduces to Y± = {Sfi)o, and Yii can also 



tain fa{z) 



~Ma zAaZ . Hence fa{z) measures the 



be calculated analytically [12|, |21|. T?he system is gapless 
below Hc2 and Y\^ , Y± are monotonic functions of H as 
shown in Fig. [TJa). 

For finite interladder couplings j > 0, the system is 
gapped for H < Hd [22| and has a. S^ = S^ ~ ground 
state. Thus we expect that Y|| and Y± show special sig- 
natures at Hd when the gap is closed and the excita- 
tions become accessible. In this range of magnetic field 
H w J±, yji can exhibit a non-monotonical behavior 
at low temperatures and strong couplings (j = 3.6, 10 



in Fig. [IJc-d)). In contrast Y± shows an intermedi- 
ate plateau like feature at low coupling (cf. j = 0.5 in 
Fig.IHb)) and a strong rise with a subsequent plateau at 
a value 0.5 for stronger coupling. 

The features at large j are roughly understood using 
a strong coupling approach. In this limit for Hd < 
H < Hc2 (gapless regime) and fc^T ^ J±, the Hilbert 
space per rung can be approximated keeping only the 
two lowest states, i.e. the triplet \t~^) = \ tt) a^nd singlet 
|s) = (I tl> - I lt>)/V2 state. Mapping these onto two 
effective spin-1/2 states |i+) -^ |t) and \s) -^ ||), the 
system H^ becomes equivalent to an effective anisotropic 
spin chain 23l426l| 



■Hxxz — J\\ /, (S'f 5'f+i + SfSf_^^ + AS-S-_^i j - HS^. 



with A = 1/2 and H = H - J± - J\\/2. Within this 
approximation 
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1/8 + {Sl)/2 + {St~St^,)/2 - jSfSf^,) 
1/2 + (5f ) 



(9) 

In particular, this approximation recovers the almost 
constant behavior of Y± in the gapless regime at low tem- 
perature. In contrast, although Y\\ « 1/2 for H ^ Hc2, 
this quantity shows a dip in the gapless regime and a 
plateau when <^ H < Hd- These effects are present 
in Fig. [IIc,d) in both the strong coupling approxima- 
tion and the full ladder at low temperature. The use of 
Eq. ([5]) for the discrimination between the two kinds of 
anisotropies J!, and J^ is very efficient in this range 
of parameters H k, J±^ and fcsT « J|| for which very 
distinct behavior contributions from Yn and Yi_ can be 
well separated. Note that when j decreases these contri- 
butions become less distinguishable (see Fig. [TJb)). 

Recently ESR measurements performed on BPCB |17| 
have pointed out small anisotropies from the perfect 
Heisenberg Hamiltonian. BPCB is known to be a very 
good realization of weakly coupled spin-1/2 ladders with 
remarkably low couplings J\\/kB = 3.6 K and J^L/ks = 
12.9 K {i — 3.6) and has been intensively investigated 
during the last decade with various experimental tech- 
niques y, |l7|, |l9|, |27H34| - This compound is composed 
of two differently oriented equivalent ladders that can be 
modelized with the Hamiltonian H. with a very low inter- 
ladder coupling of the order of 100 mK [J, [lOl, l29| which 
we neglect in the following. Using the theoretical predic- 
tion of the ESR shift ^ , we are able to extract an upper 
bound for the anisotropy along the legs and along the 
rungs assuming the limiting cases in which only 3'« and 
J^ are present, respectively. To do this we fit our results 
for the ESR shift ([5]) to the experimental data from [17| 
(Fig. Elja)) after substraction of the linear Zeeman com- 
ponent (Fig.IHb)) with the Lande factors gx = 2.09±0.01 
and gi = 2.27 ± 0.01 (for the ladders 1 and 2 labeled as 



in Rcf. [17[). In these experiments the magnetic field is 
tilted by 45° away from the b axis in the be* plan (a, b, 
c are the unit cell vectors of BPCB). 

Taking into account only the anisotropies along the 
legs J|| , the obtained fitting parameters are 

/||i/fcB = -1.5 ± O.IK , /||2/fcB=0.4±0.1K (10) 

for the corresponding orientation of each ladder respec- 
tively. Similarly, considering only anisotropies along the 
rungs J^ we get 

/±i/fcB = -1.3±0.1K , /i2/fcB = 0.3±0.1K. (11) 

Once extracted, these parameters determine the shift at 
all temperatures and values of the magnetic field. As 
shown in Fig. [21 both anisotropies are able to repro- 
duce very well all the experimental measurements avail- 
able 117|. The discrepancies between the ESR measure- 



ments and the prediction ([5]) lie within the uncertainties 
due to extraction of the peak location from the mea- 
sured ESR spectra. Indeed, the shift ([5]) is related to 
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the first frequency momentum of the peak 
stead of its maximum location which is not well defined 
when the peak is broadened. These discrepancies are 
clearly seen for temperatures T k, JajkB = 3.6 K in 
Fig- EIc) which shows the magnetic field location of the 
peak at w^ = 27r • 96 GHz extracted from extrema of the 
measured ESR spectra (Fig. 4 of Ref. [17|) versus the 
temperature. Let us note that the magnitude of the ex- 
tracted parameters pT|) considering only rung anisotropy 
is in agreement with the anisotropy/'_Li/fcB = —1.1 K 
and /±2/^B = 0.4 K extracted in [17| for a single rung 
model. However, taking the full ladder into account en- 
ables in principle the determination of the repartition of 
the anisotropy along the rung and the leg using suited 
experimental measures. 

Since the measurements in Ref. 



17[ are done in a re- 



gion where Y|| and Yx_ are not very distinct, the fit can- 
not uniquely determine their contributions. The optimal 
fit is shown in Fig. [Hb,c). More precise measurements 
around T « 3 K, for which Yn and Yj, behave very dis- 
tinctly (see Fig. [IJc)), could be used for a clear separa- 
tion of the effects of both anisotropies. The differences 
expected at that temperature are shown in Fig. [31 Note 
that one should make a compromise between the separa- 
tion of the two components Yn and Y^ and the occurrence 
of a well resolved ESR resonance peak. The extraction 
of the anisotropics J^ p and their corresponding orienta- 
tion Aq , 6q, , Ca would require measurements for different 
magnetic field directions. 

In summary, we proposed a general theory to deter- 
mine the frequency shift of ESR in spin-1/2 ladder com- 
pounds. By taking advantage of the wide applicable 
scope of the perturbative formula ([4]) , combined with nu- 
merical evaluation of the static correlation functions at 
arbitrary temperatures, we computed the shift in a wide 
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Figure 3. Theoretical prediction of the ESR shift Suj versus 
the magnetic field at T = 3 K for the two ladders orientations. 
The predictions are plotted in dotted, dashed and full lines 
considering only the rung, leg anisotropy with the fitting pa- 
rameters (|lUp . (|ll|l and with the best fitting parameters given 
in Fig. (2] respectively. The shaded area represents the stan- 
dard deviation of the ESR experiments in Ref. [ij from our 
theoretical prediction. Note that for the experimental config- 
uration of Ref. [13 , the anisotropies along the magnetic field 
for ladder 1 could be well separated from the ESR measure- 
ments at T = 3 K. 



Figure 2. (color online): Magnetic field dependence of (a) 
the ESR paramagnetic resonance frequency ujr, (b) the ESR 
resonance frequency shift Su. (c) Temperature dependence of 
the magnetic field at which the ESR paramagnetic resonance 
at ojr ~ 27r-96 GHz occurs. Symbols denote experimental data 
taken from [I3 and lines theoretical predictions of the best 
fitting parameters /||i/fcs = 0.1 K, f±i/kB = —1.4 K and 
/112/fc-B = 0.2 K, f±2/kB = 0.2 K if not stated otherwise. The 
results corresponding to the two differently oriented ladders 1 
and 2 are plotted in black and grey (online red), respectively. 
In (a) the linear Zeeman component of the ESR frequency 
resonance is represented by full lines (for the Lande factors 
gi and 32 ). We show in (b) the theoretical prediction of the 
ESR shift at r = 1.3 K for the ladder 1 considering only the 
leg (black dotted lines) or the rung anisotropy (black dashed 
lines) with the fitting parameters (jlOp and (|lip . respectively. 
These fits show only small deviations with the best fit in this 
range of temperature and magnetic field. 



range of temperature and magnetic field and showed how 
the ESR technique can distinguish the nature of differ- 
ent anisotropies. The predicted ESR shift is in very good 
agreement with experimental data for BPCB. Based on 
our results, we propose a new set of experiments to better 
resolve the different anisotropies. 
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